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ABSTRACT 



With Spitzer IRS we have obtained sensitive low-resolution spectroscopy from 5 to 35/im 
for six supernova remnants (SNRs) that show evidence of shocked molecular gas: Kes 69, 3C 
396, Kes 17, G346.6-0.2, G348. 5-0.0 and G349. 7+0.2. Bright, pure-rotational lines of molecular 
hydrogen are detected at the shock front in all remnants, indicative of radiative cooling from 
shocks interacting with dense clouds. We find the excitation of H2 S(0)-S(7) lines in these SNRs 
requires two non-dissociative shock components: a slow, 10 km s _1 C- shock through clumps 
of density 10 6 cm -3 , and a faster, 40-70 km s _1 C- shock through a medium of density 10 4 
cm -3 The ortho-to-para ratio for H 2 in the warm shocked gas is typically found to be much 
less than the LTE value, suggesting that these SNRs are propagating into cold quiescent clouds. 
Additionally a total of 13 atomic fine-structure transitions of Ar + , Ar ++ , Fe + , Ne + , Nc ++ , S ++ , 
and Si + are detected. The ionic emitting regions are spatially segregated from the molecular 
emitting regions within the IRS slits. The presence of ionic lines with high appearance potential 
requires the presence of much faster, dissociative shocks through a lower density medium. The 
IRS slits are sufficiently wide to include regions outside the SNR which permits emission from 
diffuse gas around the remnants to be separated from the shocked emission. We find the diffuse 
H2 gas projected outside the SNR is excited to a temperature of 100-300 K with a warm gas 
fraction of 0.5 to 15 percent along the line of sight. 



Introduction 



As a supernova remnant evolves into its surrounding environment , it becomes the dominan t reservoir of 
energy that drives dynamics and chemistry in the interstellar medium ( McKee &: Ostrik"er1ll977 ). As massive 



stars in general do not drift far from their parent cloud in their short lifetimes, supernovae near giant 
molecular clouds (MCs) should be common. Yet there are only a handful of SNRs which have been clearly 
identified as interacting with MCs. Pure-rotational transitions of molecular hydrogen(H2) are excellent 
indicators of such interactions as high columns of H2 can survive the passage of the shock. 

Molecular hydrogen is the most abundant molecule in the Universe, playing a crucial role in the interstel- 
lar mcdium(ISM) as a catalyst of chemistry and a major coolant. H2 is excited by energetic processes such as 
shocks and photo-dissociation regions, with numerous bright emission lines in the near- and mid-infrared (IR) 
available as detailed diagnostics of the excited gas. The low-lying pure rotational transitions of H2 probe 
the bulk of the gas involved in the shock interaction and are excellent diagnostics of the shock conditions. 
Due to their relatively low critical densities the S(0) to S(7) transitions are not good diagnostics of the ex- 
citation mechanism; however, supernova remnant shocks have the advantage of being comparatively simple 
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systems in which to understand detailed shock processes as there is no complicating stellar radiation source. 
The study of Galactic supernova remnants serves as an important temp late for other shock in teractions, for 
example the spectacular galaxy-scale shocks seen in Stephan's Quintet ( Appleton et al. 20061 ). 



In this paper we report the results of mid-IR long-slit spectroscopy of six SNRs which are identified as 
bright in IR molecular lines. These SNRs have not been previously studied in the mid-IR and represent a 
diversification of the common sample of nearby ineracting SNRs such as IC 443, W28, W44 (see Rho et al. 
2001 and Reach & Rho 2005 and references therein). The observations are briefly described in §2, followed 
by a summary of the properties and observed line emission for each source in §3. Section 13.31 pertains to 
the physical conditions of the detected pure rotational H2 lines. Sections 14. 1 1 and [ 4.21 reconciles the observed 
molecular and atomic lines, respectively Section [5] presents a discussion of diffuse H2 observed outside the 
shock front. 



Observations 



Motivated by the detection of eighteen IR-bright SNRs using the GLIMPSE survey |Reach et al.ll2006[) . 



we have selected six SNRs which display Spitzer Infrared Array Camera (IRAC) colors indicative of molecular 
shocks. Of the six interacting remnants we consider in this work, four SNRs - Kes 69, G346.6-0.2, G348. 5-0.0, 
and G349.7+0.2 - have previously identified 11(1720 MHz) masers wh ich have been established as excellent 
signposts of interaction with an adjacent MC jYusef-Zadeh et al.ll2003h . For Kes 17 and 3C 396 the evidence 



for interaction was l ess clear, but both have bright filamentary radio shells and IRAC band ratios indicative 
of molecular shocks ( Reach et al. 20061) . 



The Spitzer Space Telescope Infrared Spectrograph (IRS) was used to obtain spectra of these SNRs to 
characterize the nature of emission using the short-low(SL, 5.2-15/im) and long-low(LL, 14-42. 4/im) modules 
and with a resolving power of R « 60-125. In each remnant a single position coincident with a peak in the 
IRAC channel 2 emission was chosen (see Table [1]). Exposure times of 150 and 180 seconds were used for the 
SL and LL modes. As the exposures were nodded, the effective integration time for the overlapping region 
between both orders is doubled typically producing a line signal-to-noise ratio of 20. For G349.7+0.2 the 
remnant is bright enough that only half the integration time was needed. A cluster of reference sky positions 
devoid of IR sources were observed within a few degrees of the targets. Subtracting the median averaged 
reference position from the source positions successfully mitigated most rogue pixels. The remaining rogue 
pixels were removed manually utilizing the IRS CLEAN program provided by the Spitzer Science Center. 
Additionally the slits were nodded with respect to the target position so that different sections of the slit 
observed the same overlapping target position. 

Spectra were extracted from the region of overlap between the IRS SL and LL slits and aperture 
correction was applied. We did not find it not necessary to defringe the observed spectra. For each target 
we subtracted a local background using a nearby position within the same IRS nodded pair. In all cases a 
position which lies outside the remnant, and has the lowest background region from one of four co-added 
observations is used. As the short-low and long-low slits are aligned roughly perpendicular to each other, we 
must use different locations as representative local backgrounds for the short- and long-orders. In practice 
the backgrounds did not differ substantially and produced only minor discontinuities between the two IRS 
orders. The H2 S(2) through S(7) lines are only detected from the SNR interaction regions and are not 
affected by a local background subtraction. However, the H2 S(0) and S(l) lines show diffuse emission across 
the slit, even outside the SNR, but with a clear enhancement of emission at the SNR shell. The detailed 
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properties of the remnant and diffuse background are discussed separately in ^3.31 and £j5j respectively. 



3. Results 

3.1. IRS Spectra and Line Brightnesses 

IRS spectra of SNRs Kes 69, 3C 396, Kes 17, G346.6-0.2, G348. 5-0.0 and G349. 7+0.2 after local back- 
ground subtraction are presented in Figure [JJ with prominent emission lines and PAH features labeled. For 
all remnants bright pure rotational transitions of H2 (0,0) S(0) through S(7) are detected. Atomic sulfur and 
several ionized fine-structure transitions - Ar + , Ar ++ , Fc + , Ne + , Nc ++ , S ++ , and Si + - are also observed. 
The multiple molecular and ionic transitions detected with Spitzer are the first clear evidence for interac- 
tion in most of these remnants. For 3C 396 there is evidence that molecular and ionic lines are spatially 
segregated within the IRS slits shown in Figures [2] and [3l A steadily rising background continuum level 
and in some cases PAH emission are also seen after local background subtraction, but detailed analysis and 
discussion of this emission is deferred to a future paper. In general there is good agreement between the 
IRS spectra and the nature of IR emission predicted from the simple color template analysis of IRAC band 



emission by iReach et al.l (|2006l ) 



Properties of the SNRs considered in this work are given in Table [TJ The best distance determination is 
used to obtain a corresponding physical size of the remnants radio shell. For Kes 69, 3C396 and G349.2+0.2, 
we used extinction values from the literature. For Kes 17, G348. 5-0.0 and G346.6-0.2, we computed the 
visual extinction in two or three methods listed below and found the most consistent values. First, the IRS 
spectra are fit by the lines and continuum (using PAHFIT; Smith et al. 2007) where the silicate absorption 
dips are important to determine the extinction. We have verified the extinction solutions of PAHFIT by 
estimating the de pth of the 9.7 /xm silicate feature against the continuum of the remnant. We used the 
extinction curve of Draine ( 20031 ) with a visible extinction parameter, Ry = Ay/E(B-V), of 3.1 appropriate 



for the diffuse ISM in the Milky Way to derive the extinction from the 9.7/im optical depth. Second, the 
total hydrogen column density along the line of sight is derived from X-ray observations for Kes 17 and 
G348. 5-0.0 using archival ASCA data. The detailed X-ray analysis will be presented in subsequent paper 



(Rho et al. in preparation). Third, we have used CO and HI data (jDickev fc Lockmanlll990l : I Dame et al 



20011 ) to estimate the total hydrogen column to the SNR. When the kinematics of the SNR are known from 
masers, we can estimate the column of atomic and molecular hydrogen by integrating up to the systemic 
velocity of the SNR. This gives an upper limit to the hydrogen column column which is appropriate if the 
SNR is located at the far distance solution of the Galactic rotation curve. This was an important constraint 
for G348. 5-0.0 as it is the only SNR without detected central X-ray emission or an associated pulsar, and 
it is believed to lie at the far distance near SNR CTB 37A. Given there is still some uncertainty in the 
appropriate de-reddening values for these SNRS we list the observed line intensities after local background 
subtraction in Table 2; the values after extinction correction are given in Table 3. 

Using the distance and extinction given in Table [TJ the total luminosity of IR lines observed for Kes 69, 
3C 396, Kes 17, G346. 6-0.2, G348. 5-0.0 and G349. 7+0.2 are derived. In all cases the lowest luminosity source 
was Kes 69 and the brightest was G349.7+0.2. For H 2 S(0) to S(7) we find a total luminosity of 1.4 to 68 L Q . 
For ionic lines with ionization potentials less than 13.6 eV (Fe + , Si + ) and greater than 13.6 eV (Ar + , Ar ++ , 
Ne + , Ne ++ and S ++ ) we find a total luminosity between 5 to 1212 and 0.8 to 595 Lq respectively. These 
two groups of ions arc considered separately because they require different physical condition s and previous 



spectral line mapping of interacting SNRs has shown the two are not spatially correlated (jNeufeld et al 
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20071 ). Line luminosities for each species are given in Table [4] for comparison. The [Si II] 34.8/xm line is 



always observed as a significant cooling line. For Kes 69, Kes 17 and G346.6-0.2 molecular hydrogen lines are 
also significant coolants in the observed IRS wavelengths. For 3C396 and G348. 5-0.0 the cooling from the 
three main ionic lines of Fe, Ne and Si are comparable to the cooling from H2 lines. The youngest remnant, 
G349.7+0.2, has bright ionic lines that are far more luminous than H2. We note that the luminosities given 
are for the observed apertures which varies with wavelength. For the IRS SL and LL modules the apertures 
are 3.6 x 8"and 10.5x20" respectively. 



3.2. Individual Supernova Remnants 



Here a brief description of the properties of each SNR is given, noting the existing evidence for interaction 
and the individual results of our Spitzer IRS observations. 

G21.8-0.2 (Kes 69): The SNR Kes 69 has a clear Southern radio shell but little radio emission to the 

North. IRAC 4.5^m images, that include emission from H2 show a faint filament along the prominent 

Southern radio shell. A single OH(1720 MHz) maser with a velocity of +69.1 km s _1 has been detected in 

Recent CO and HI observations show Kes 69 is 
1 



the far Northern extent of the remnant (jFrail et al.lll996l ) 



at a distance of 5.2 kpc with a systemic velocity of +80 km s 1 (jTian fe Leahvll2008r). Near th is velocity 
OH(1720 MHz) emission along the radio shell indi cates a molecular shock |Hewitt et al" 2008h . ROSAT 



observations give a column of ATh =2.4x10 cm ([Yusef-Zadeh et al 



20031) . 



G39.2-0.3 (3C 396): 3C 396 has a bright radio shell and a central X-ray emission. The western shell 
appears brighter in both radio and X-ray emission, suggestive of interaction with a MC. IRAC images in all 
bands show a bright IR filament along the western edge with colors indicating an ionic or molecular shock. 
There is a compact centra l X-ray nebula with a photon index of 1.5 consistent with a central pulsar wind 



nebula ( Olbert et al. 2003h . There is also significant therm al X-ray emission from the interior of the SNR 



not related to the pulsar wind nebula ( Harrus fe Slanel 1999h . The total hydrogen column for the supernov a 



remnant measured by both Chandra and ASCA is 4.7x 10 22 cm 2 |Harrus fe Slane!ll999tlo"lbert et al.ll2003h 



Assuming the SNR and pulsar wind nebula are related, a distance of 6.5 kpc is obtained from the hydrogen 
column. 

The obtained spectrum of 3C 396 from the western IR filament, shows a wealth of both H2 and ionic 
lines of Ar, Fe, Ne, S, and Si. Figure [2] shows a clear spatial separation between the H2 and ionic lines along 
the western shell. This is clear evidence for spatial variation of physical conditions along the shock front. In 
Figure [3] diffuse interstellar emission lines of H2 S(0), [S III] and [Si II] are seen outside the SNR, whereas 
[Fe II] is only prominent inside the remnant. These long slit observations hold an advantage over spectral 
mapping of nearby remnants in that the emission associated with the SNR can be clearly distinguished from 
diffuse emission using local background subtraction. 

G304. 6+0.1 (Kes 17): An inco mplete radio shell, K es 17 has prominent IR filaments with colors indicating 



suggests a distance of at least 9.7 kpc (jCaswell et al 



bright emission from shocked H2 (IReach et al. 2006). Th e remnant is not well studied, but HI absorption 



1975). Given the absence of clear kinematic signatures 



we give lower limits for the distance and size of Kes 17 in Tabic [TJ Archival ASCA observations show diffuse 
interior X-ray emission, allowing us to derive a total hydrogen column of 3.6xl0 22 cm~ 2 toward Kes 17. 

G346.6-0.2: This remnant has a faint radio shell filled with diffuse thermal X-ray emission, making it a 
mixed- morphology SNR. Five OH(1720 MHz) masers are detected along the Southern shell of the SNR near 
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a vel ocity of -76 km s" 1 ; Z eeman splitting of the maser line gives a line-of-sight magnetic field of +1.7±0.1 
mG ( Koraleskv et al. 19981) . The kinematics implied by the velocity of the associated maser gives a near/far 
distance 5.5 or 11 kpc. Archival ASCA GIS observations show diffuse X-ray emission from the remnants 



interior which is best fit by a total hydrogen column of 2.68^Q gg x 10 22 cm 2 and a temperature of kT 



1.36^o'36 e V- The large hydrogen column and small size of the remnant favor the further distance of 11 kpc, 
which we adopt here. We note that MIPS 24^m observations show diffuse emission covers the remnant, 
likely due to a foreground HII region to the Northwest. 

G348. 5-0.0: An overlapping radio shell detected adjacent to CTB 37A, the SNR G348. 5-0.0 appears to 
be a kincmatically distinct interacting remnant. Both CO and OH maser observations in dicate interaction 
with a cloud at -20 km s" 1 for G348. 5-0.0, whereas CTB37A is found at -65 km s" 1 jFrail et all 



Revnoso fc Ma 



ngum 



19961 : 



2000). A line -of-sight magnetic field strength of -0.8±0.1 mG is measured for a maser 



in G348. 5-0.0 (|Brogan et alJl2000h . The near/far distance ambiguity for both G348.5-0.0 and CTB 37A are 
resolved in favor of the far distance by the HI observations, but not with sufficient resolution to differentiate 
between the two bright radio sources (jCaswell et al.lll975l ). Archival ASCA observations did not allow a 
determination of the hydrogen column density towards G348. 5-0.0 so using the established kinematic velocity 
of -22.0 km s^ 1 we integrated the total CO and HI column densities to estimate a column of 2.8 xlO 22 cm~ 2 
towards the remnant. 

G349.7+0.2: The most luminous and distant SNR within the galaxy with detected OH(1720 MHz) masers 
is G349.7+0.2. Numerous molecular transitions indicative of shock interaction are detected toward includ- 
13 CO, 12 CO, CS, H 2 , HCO+, HCN, H 2 CO and SO |Lazendic et alJl2004h . A kinematic distance to 



mg 

G349.7+0.2 places it at 22 kpc, but rotation curves at such a large Galacto-centric radius are not well con- 
strained. X- ray observations yield a neutral hydroge n column of 5.8xl0 22 cm -2 and suggest an age of only 



3,500 years (jLazendic et al.ll2005l) . iFrail et al.l (|1996l ) identified five maser spots at velocities between +14.3 
and +16.9 km s _1 . The impacted cloud appears at fT,s R=+16.2 km s" 1 wi th a linear size of about 7 pc, a 
mass of ~104 M Q and a volume density of ^10 3 cm -3 (jDubner et al.ll2004l ). The SNR-cloud interaction is 
thought to be taking place on the far side of the SNR from a comparison between total absorbing column 
density and slightly red-shifted wings in the optically thin 13 CO spectra. Figured] shows the IRS LL2 spectra 
across the remnant. Bright H 2 S(l) and [Fe II] lines are seen only associated with the SNR shell. Though 
some [Ne III] and [S III] is present from diffuse interstellar emission outside the remnant, both lines arc 
clearly enhanced coincident with the bright shell. 



3.3. Excitation of Molecular Hydrogen 

The bright emission from the H 2 S(0)-S(7) lines present in all six SNRs is an excellent diagnostic of the 
physical conditions of the shocked gas, of which H 2 is the dominant constituent. The energy level diagram 
for H 2 is presented in Figure [5] with the transitions observed by Spitzer labeled (by conventional branch 
notation S indicates A J = -2). As these lines are optically thin they provide a direct measure of the column 
density for each upper state J u = 2 to 9. We note the S(7) transition has a rotational energy higher than 
that of the bright H 2 1-0 S(l) line typically observed in the NIR. 

Rotational diagrams from dereddened H 2 line intensities are presented in Figure [6l Here the column 
density per magnetic substate is plotted against the upper energy level, such that a single temperature 
distribution in LTE follows a straight line. It is clear that a simple LTE excitation model for the observed 
level populations cannot fit both the positive curvature resulting from a range of gas temperatures, or the 
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zig-zag pattern between subsequent ortho- and para-states out of equilibrium. However, it is instructive to 
examine the properties of a two temperature model which can be well fit to the data. Table O gives the 
best fits from least squares fitting assuming two temperature distributions and simple LTE excitation of H2 . 
Common for the fits are a warm component with T w ~ 250-600 K and a column N w ~ 10 20 cm -2 , and a 
hot component with ~ 1000-2000 K and a column N„ ~ 10 19 cm -2 . In Figure[6]the warm component is 
drawn as a dashed line, the hotter component is drawn as a dotted line and the total contribution as a solid 
line. For the hotter gas component (dotted line) the ortho to para ratio (OPR) is found to be at or near the 
statistical equilibrium value of three. However, the warm component (dashed line) shows a range of OPR 
between 0.4 and 3, manifested as zig-zags between ortho- and para- levels in the energy diagrams. This larg e 
range in OPR™ was also found for the interacting SNRs 3C391, IC443, W28 and W44 jNeufeld et alj|2007l) . 
The conversion of ortho to para in shocks is strongly dependent on temperature, w ith threshold temper atures 
at 700 K to start conversion and 1300 K to rapidly reach the equilibrium value (jTimmermannl 119981 ). The 
temperature of warm gas is substantially less than 700 K for all but SNR G346.6-0.2 which has an OPR lu 
of 3. For the other remnants an OPR between 0.4 and 1.4 is characteristic of an equilibrium temperature 
between 50 and 100 K. For shock velocities < 20 km s -1 the H2 gas w ill retain the OPR of the pre-shock gas, 
with only a small dependence on gas density (jWilgenbus et al.ll2000| ). Low OPR values indicate that there 
has not been substantial prc-heating of the pre-shock material with which these SNRs are now interacting. 



4. Shock Models 

At the position of the bright IR peak observed in each SNR both strong molecular and ionized species 
are present. This discordant mix of spectral lines cannot be reconciled with a single shock into a uni- 
form medium. We consider two general classes of shock models. Jump(J) shocks have a discontinuous 
change of hydrodyanmic var iables, with heating over a negligible thickness sufficient to dissociate molecules 
(|Hollenbach fc McKed[l989l) . [t is possible for H2 and other molecules to reform in the cooling post-shock 



region. In contrast continuous(C) shocks have a transition region in which friction between ions and neutrals 
occurs over a sufficiently large thickness that the t ransition from pre- to post-shock is continuous and nondis- 
sociative for molecular gas of modest ionization ( Draine et al. 19831 ). In the following sections we consider 



the physical parameters of shocks which can explain the observed IR line emission. 



4.1. Ho Rotational Lines 



For each remnant the observed H2 columns are compared to published C- and J-shock models. The H2 
line intensities, particularly at low excitation temperatures, are typically more than an order of magnitude 
brighter than predicted by J-shock models (eg., Hollenbach & McKee 1989). A comparison with C-shock 
models (eg., Draine et al. 1983) shows lines of comparable magnitude, but a single C-shock could not 
simultaneously account for the strengths of all observed H2 lines. 

Using a grid of shock models a detailed fit to the excitation of H2 lines was performed using least squares 
fitting. The relevant shock variable s are shock veloc i ty v, s and pre-shock density n D . Two recent C-shock 
models were used. The first model of lWilgenbus et al.l (|2000T ) accounts for the non-equilibrium OPR observed 
in the warm H 2 gas. A grid of computed models spans (logio n D )=3,4,5,6 cm~ 3 , v s =10,15,20,25,30,40 km 
s _1 and OPR=0.01,1,2,3. For shock velocities above a certain critical velocity H 2 is dissociated; without H 2 
cooling a continuous transition cannot be sustained and the shock becomes J-type. iLe Bourlot et al.l (|2002l ) 
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found that C-shocks could be sustained at much higher velocities than previously t hought, ranging up to 2 3 
to 80 km s _1 for densities of 10 7 to 10 3 cm" 3 , with i ncrements matching those of IWilgenbus et al.l (|2000f ) . 



The second shock model from lLe Bourlot et al.l (|2002i ) correctly accounts for shock speeds up to the cirtical 
velocity but assumes OPR=3. However, for shock v elocities much less than the critical velocity (~ 10-20 
km s" 1 ) there is good agreement between the work of Le Bourlot et al. ( 2002h and Wilgenbus et al. ( 200fjf ) 
for OPR=3. Together these two sets of models form a grid of C-shock models that span the observed H2 
excitation parameters found in Section [3.31 



In our model fitting each shock component is given a free linear scaling to account for uncertainty in 
the geometry and area of the shock relative to the IRS aperture. This factor <E> S is the ratio of the projected 
beam size to the surface area of the shock. If the beam is not completely covered by shocked material the 
value of <fr s can be less than unity. 

No single shock model is found that can reproduce the positive curvature observed in the Boltzmann 
diagrams. For a good fit to the data a combination of two shock models is required. A slow C-shock gives a 
good fit to the lower excitation lines (eg., H 2 S(l)) but under-predicted the observed b rightness of higher ex- 
citatio n lines (eg., H 2 S(7)). In addition to C-shocks we included the J-shock models of Hollenbach fc McKeel 
(|l989l ) as a discontinuous grid of models with densities of 10 3 -10 6 cm -3 and shock velocities of 30-150 km 
s" 1 in 10 km s" 1 increments. The results using a single C-shock with either a second faster C-shock or a 
fast J-shock are given in Tables [6] and [7l Figure [7] shows the best fits to the observed excitation of each SNR 
which are generally found to be two C-shocks. 

We find a slow ~10 km s" 1 , dense ^10 6 cm" 3 C-shock gives an excellent fit to the S(0)-S(3) lines, 
including non-LTE OPR. The higher excitation lines arc well fit by a fast C- or J-shock through a lower 
density medium (n Q ~10 4 cm -3 , v s ~30-60 km s" 1 ) in which the OPR is quickly brought to equilibrium. 
In general, a fast C-shock gives a better fit to the observed excitation of H2 than a fast J-shock. The 
only exception is G346.6-0.2, in which a single C-shock at a velocity of ~15 km s _1 into a density of 10 5-6 
cm" 3 and an OPR of 3 reproduces the observed H 2 S(0)-S(6) lines but underestimates the H 2 S(7) line. 
Either a fast C-shock or a fast J-shock into a lower density medium can produce the S(7) line but this 
shock component is not well constrained so either solution is viable. Figure [8] shows the alternative J-shock 
solution given in Table [8] For all other remnants a C- and J-shock combination do not do a good job or 
reproducing the observed excitation, however we do list the best-fit parameters in Table 7 for comparison 
with the two C-shock models. It is p ossible that our best fit of two C-shock models is a characteristic of 
temporal evolution of an MHD shock ( Chieze et al. 1998t Cesarskv et al. 1999 ). 



We compare the shock fits to a simple analytic treatment of energy conservation in C-shocks (jNeufeld et al 



2006L Appendix B) in order to check the validity of fitting results and to understand associated physical pro- 
cesses. The characteristic gas temperature of a C-type shock is given by 



T s = 375 b" a3b (v s /10 km s 



-IU.35 



K. 



(1) 



where b = B/fiG (n /cm" 3 )" - 5 . Given the standard assumption that b=l the warm and hot gas components 
discussed in Section [3~3l typically correspond to shock velocities of ~10 and ~30 km s" 1 , respectively. This is 
in general agreement with the parameters fit in the modeling, however the hot gas shock velocity is typically 
found to be somewhat higher. The observed column density of shocked H 2 is given by 



4 x 10^ b n^ 4 5 (vs/10 km s) 



-0.75 



4's 



(2) 



where the preshock density is 10 4 n 4 cm 3 and <i> s is the ratio of the projected beam size to the surface area 
of the shock as before. Even though the preshock density, n Q is not well constrained by the low-J rotational 



- 8 - 



transitions, a comparison between the H2 column for the warm and hot component in Table [5] and the 
expected value of N s /$ s from Equation[2]can establish whether the fitted <£> s is reasonable. Using the best-fit 
model preshock density and shock velocity in Table [5] we find the expected values of $ s match the model- 
fitted $ s within a factor of a few. Typically $ s is found to be 0.1-0.5 and 0.03-0.02 for the warm-slow and 
hot-fast shocks, respectively. For Herbig-Haro objects there is significant substructur e not resolved by IRS 



and values of <&. s for the w arm and hot components are typically 0.2ng 4 5 and 0.05n,Q4 5 (jNoriega-Crespo et al 



2002 



Neufeld et al.1 120061 ). Supernova remnants also appear to have significant substructure which IRS 



observations are not able to resolve. 

We now compare the shock models of that fit of our H2 lines with other models and observations. Le 
Bourlot et al (2002) applied C-shock models to the H2 lines from Orion and found a best fit with n o =10 4 cm -3 
and two different shock velocities of 40 and 60 km s~ 2 . Table M shows that two different density models yield 
a better fits than two different velocities for a single density. iBurtonl (J1992J) suggested a density of 10 6 cm -3 
with a very low velocity of 10 km s _1 ; such a model cann o t expl ain all the higher- J observed H2 emission (see 
the dashed and dotted lines in Figure [7]). Richter et al. ( 19951) suggests partially dissociated H 2 , from a J- 
shock, but the implied high velocities could not be explained. Fits to ISOCAM ob servations of IC 443 s how a 
density of 10 4 cm~ 3 with excitation temperatures of 330-657 K and 1115-1300 K (jCesarskv et al.lll9 99). Our 
fit results show two main differences from those of Cesarsky et al. (1999); one is a higher density component 
and the other is departure from LTE OPRs. The main differences may arise from our inclusion of S(0) and 
S(l) lines which could not be observed by ISO. The departure from LTE is independent of shock models and 
implied by the two temperature fitting of H2 excitation diagram. The inferred OPRs are out of equilibrium 
i ndicating the time-scale of the sh ock passage was too short to establish an equilibrium ortho/para ratio 
(jPineau Pes Forets fc Flowerll2000f ) . In the course of shock passage, temporal non-equilibrium could be due 
to dynamical processes with colder gas, or fast conversion of ortho-H2 to para-H2 during the formation on 
the surface of grains. Molecular hydrogen formation is mainly governed by the properties of grain surfaces 
and shocked medium quickly accreted in t he diffuse ISM in even a modest velocity (^30 km s" 1 ) shock 
(jjones et al.1 Il994t ICazaux fc Tielend 120021 ) . A combination of two C-shock models indicates a shock into 



very dense (10 cm 3 ) gas with a low (10-15 km s ) shock velocity where sputtering and dust destruction 
may not be as efficient. 

Given a multi-phase ISM a combination of pre-shock densities is possible. A modera te density (10 3 -10 4 



cm 3 ), faster shock component has been resolved in W44 and W28 (jR.each fc Rhol l 20001 ) . Widths of H2 m 
the shock front are measured to be ~10 16 cm, consistent with nondissociative shocks into gas with densities 
of ^10 4 cm -3 . These densities are consistent with our fast C-shock component of 40-70 km s _1 (see Tablc[B]). 
Perhaps surprisingly the slower, denser shock has a larger filling factor than this faster, lower density shock in 
our fitting. At such high densities of 10 5 -10 6 cm -3 cooling will be faster with a lower peak temperature and 
a larger shock thickness if the two shocks are roughly isobaric. Cooling could be significant from molecular 
lines other than H 2 (such as H 2 0) though the current observations are not sufficient to detect such lines, for 
which observations with a higher resolution spectroscopy arc required. 



4.2. Ionic Lines 

Strong ionic lines observed toward all SNRs and cannot be reproduced by steady-state C-shock models 
(see Draine et al. 1983). We use the J-shock models to derive the shock velocity and pre-shock density 
purely from consideration of the surface brightness of observed ionic lines. Here our purpose is to determine 
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the gross properties of the fast J-shock responsible for the observed fine-structure linesQ. 

Diagnostics of the [Fe II] lines use excitation rate equations which are presented in Rho et al. (2001). 
Figure [5] shows contours of line ratios of 17.9/5.35/im and 17.9/26/xm; the former ratio is mainly sensitive 
to density and the latter ratio changes depending on both density and temperature. Updated atomic data 



(jRamsbottom et al.l 120071 ) affected (at the <10% level) the contours at low densities (<200 cm -3 ). The 
estimated density and temperatures from [Fe II] ratios are summarized in Tabled A set of electron density 
and temperature solutions are possible for each of G349.7+0.2 and 3C 396: (n e , T) = (700 cm -3 , 7000 K) 
and (270, 2.3xl0 4 ), respectively. The ratios of Kcs 17 indicates (n e , T) = (100, 7000) showing a low electron 
density. For Kes 69 and G348. 5-0.0, the two ratios did not converge; the ratio of [Fe II] 17.9/5.35yum indicates 
a higher temperature than that from the ratio of [Fe II] 17.9/26/im. This may be due to either uncertainties 
of collisional strengths at low density and temperature or may indicate presence of multi-temperature gas. 
The ratio of [Fe II] 17.9/5.35 /jm indicates a low density (< 400 cm -3 ) and low (<2000 K) temperature. 
Temperature could not be well-constrained; lines of higher upper energy levels such as the [Fe II] 1.64/im 
line are required to constrain temperature. 

The pre-shock medium density n D and shock velocity v s are also constrained using the strength of the 
[Ne II] 12.8, [Si II] 34.8, [Fe II] 26/im lines (called the three main ionic lines hereafter, see Figure 6 of Rho et 
al. 2001 which were generated based on Hartigan et al. 1987 and McKce et al. 1984). The inferred medium 
densities and shock velocities are summarized in Table [9] Iron lines have a high ionization potential and 
may not arise from the same physical conditions as these ionic lines with lower ionization potentials. 

The ratio of [Ne III] 15.5/[Ne II] 12.8 is also useful to determine shock velocity and is presented in 
Figure [TO] with results summarized in Table [5] The inferred shock velocities from Ne ratios are generally 
consistent with those from the brightnesses of three main lines except 3C 396. The only remnant for which 
the Ne III/Nc II ratio is greater than unity is 3C 396, suggesting a shock in excess of 300 km/s while the 
inferred shock velocity from the line intensities of [Si II], [Fe II], [Ne II] are 100-150 km s _1 . The [Ne III] in 
3C 396 may arise from fast shocks than the lower-ionization state lines. 



4.3. Molecular and Ionic Shocks and Ejecta Contribution 



Estimates of the shock ram pressure (P s ) can be derived from the shock modeling in the previous two 
sections. The best-fit parameters from both C-shocks derived from H2 and the fast and diffuse J-shock traced 
by ionic lines are given in Tables [6] and [Qj respectively. The highest pressures are generally experienced in the 
molecular C-shocks in which the post-shock gas cools and condenses and can obtai n a pressure higher tha n 
the thermal pressure of the remnant proportional to the gas compression factor (jMoorhouse et al.lll99ll ). 
It is well established theoretically that pressure e nhancements of 10-100 times can occur when the SNR 



shell interacts with dense clumps (| Chevalierl 1 1 999^ . OH maser emission detected from Kes 69, G348. 5-0.0, 



G346.6-0.2, and G349.7+0.2, which is evidence of the survival of higher density clumps. We caution that the 
pressures derived from shock models may have uncertainties up to an order of magnitude. Given this large 



1 The [Fe II] 26/jm is prominently detected but the spectral resolution of IRS is not sufficient to distinguish between the [O 
IV] 25.9 and [Fe II] 26/im lines. Observations of interacting remnants with ISO and IRS at high spectral resolution are able to 
separate the two lin e s, and find that the [O IV] 25.9/xm line does not contribute significantly if at all to [Fe II] 26/xm line flux 
teeach fe Rholll999l ; iNeufeld et alj|2007t ). Furthermore the [O IV] 25.89/xm line has an ionization potential of 54.9 eV and is 
only detectable for very fast shocks in excess of those observed here. We expect that the [Fe II] 25.98^tm line is by far dominant 
for most if not all of the SNRs we consider here because other transitions of [Fe II] lines are detected. 
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uncertainty it is still clear that H2 emission originates from regions where the pressure is of order 10 _7 ~10~ 6 
dyn cm -2 . The pressures derived for each of the two C-shock components is consistent with an isobaric 
shock into a multi-phase cloud with a range of densities. In general the pressure derived from the three ionic 
lines are comparable to that derived for the C-shocks, though the densities are lower (100-1000 cm -3 ) and 
shock velocities higher. Again this is consistent with the SN R shock interacting w ith a multi-phase MC with 
both an ambient molecular phase and a dense clump phase ( Reach fc Rho 2000f) . 



There are a few exceptions to the general observations above. G349.7+0.2 has an exceptionally high 
ionic shock ram pressure, even exceeding that of the dense C-shocks. The bright ionic lines, inferred high 
shock velocity and a high electron density inferred from [Fe II] line suggest that this SNR is younger than 
other SNRs in our sample, and may contain a contribution from SN e.jecta. X-ray ob servations also show 



abundance enhancements indicative of a significant ejecta mass (jLazendic et al.l 120051 ) and the esatimated 
age of the SNR is 3500 years. 

Another exception is G346.6-0.2, where only very weak ionic lines and OPR of 3 are observed. The 
ionic pressure is estimated to be only ~10~ 8 dyn enr 2 whereas the pressure from the C-shocks are both 
4xl0 -6 dyn cm~ 2 . Given that the OPR is observed at equilibrium, this may suggest that the shock has had 
sufficient time to establish a steady-state C-shock, whereas those remnants with non-equilibrium OPR and 
observed ionic emission may have not yet obtained steady-state. The weak ionic lines appear to arise from 
a shock into a low density ambient medium of 100 cm~ 3 and not the dense clump material from which H2 
emission originates. 

Finally, 3C 396 shows a large range of pressures. The pressure of the hot thermal X-ray ga s observed 
in the interior of the SNR is estimated to be of order 3xl0 -9 dyn cm -2 ( Harms fe Slanel 1999 ). This is a 



factor of a few smaller than the shock pressure derived from ionic line diagnostics 1.8xl0 -7 dyn cm~ 2 . In 
comparison to the observed H2 emission from 3C396, the shock pressure is observed to be a factor of ^10 
higher. Given the spatial segregation of the ionic and H2 lines, it is plausible that H2 emission is arising 
from interaction with dense clumps while emission from the main ionic lines arises from the same shock 
propagating into a lower ambient density medium or contribution from ejecta. The fact that Fe line emitting 
gas is possibly hotter and denser than the other four SNRs (see Table[9|) and the two different shock velocities 
support that the infrared emission from 3C 396 is composed of both ejecta fragments and shocked ISM. Full 
spectral mapping (such as of Cas A; Rho et al. 2008) would be required to separate distribution of ejetcta 
and shocked ISM. 



Diffuse Warm Hg 



Warm diffuse H2 gas is known to be abundant throughout the galaxy (|Grv et al.ll2002t : iFalgarone et al 



20051 ). Diffuse H2 emission has also been noted at t he SNR-cloud interaction site, though it's relation to 
the shock interaction is unclear (jNeufeld et al.l 120071 ) . Using careful local background subtraction we are 
able to isolate the H2 emission due only to the supernova remnant interaction. We find enhanced S(0) and 
S(l) emission from the shock interaction, but often the diffuse emission is dominant. Here we discuss the 
excitation of this warm diffuse molecular gas. 

Diffuse S(0) and S(l) lines are detected for nearly every sight line, including our sky reference positions. 
The diffuse component we observe could arise from conditions in the pre-shock clouds or from low density 
diffuse gas present along the line of sight through the Galaxy. Several excitation mechanisms for warm diffuse 
H2 present in the Galactic plane have been proposed including the Galactic radiation field, H2 formation 
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energy (Spitzer 1974), and MHD shocks or small-scale turbulence (Falgarone et al. 2005). However, higher 
energy transitions such as S(2) through S(7) are not detected as for the SNRs. Because we have only one 
ortho and one para line, the ortho-to-para ratio and gas temperature are degenerate so we can only constrain 
the excitation of diffuse H2. 

Table [10] gives the fitted parameters to the warm diffuse H2 emission in the regions used as a local 
background in data reductions. We have not included Kes 69 or Kes 17 as the slit of one of the LL orders 
lies almost entirely inside the SNR. As the temperature is degenerate with the ortho-to-para ratio, we have 
assumed a range of different OPR. In the first block of physical parameters listed in Table 10 we assume 
an OPR of 3.0 to derive excitation temperatures that range from 100 to 160 K, with the notable exception 
of G346.6-0.2. Figure [Til shows the Boltzmann diagram for diffuse emission from outside the SNRs. The 
slope of the line connecting the S(0) and S(l) transitions is markedly different for G346.6-0.2 than for the 
other remnants. Diffuse emission towards G346.6-0.2 may be associated with a foreground Hll region which 
appears to dominate the MIPS 24^m image of the SNR's vicinity. 

The assumption of an ortho-to-para ratio of three gives the minimum temperature of the gas (see Table 
10, row 1). At such low temperatures the equilibrium OPR is roughly two and may be a more appropriate 
value. However, if the diffuse warm H 2 has only been recently heated from a colder, quiescent state it could 
have a significantly lower OPR. For comparison with the shocked gas, we use the OPR of the warm shocked 
gas listed in Table O These solutions for various ortho-to-para ratios are all given in the second block of 
parameters listed in Table ITDl 

Even for the lowest OPR (the third block of parameters listed in Table 10), which is taken to be the 
same as that observed for the SNR shock, the excitation temperature of the low-lying rotation transitions 
is comparable to the 180-390 K excitation temperatures observed by Savage et al. (1977) in a study of UV 
absorption sight lines. Falgarone et al. (2005) reports on the excitation of S(0) through S(3) transitions for 
a line of sight devoid of massive star forming regions finding comparable surface brightness for the S(l) and 
S(0) transitions. As we did not detect the S(2) through S(7) lines for the diffuse molecular gas we were not 
able to probe these proposed signatures of dissipative turbulence. 

Though the excitation conditions are somewhat uncertain, we are able t o estimate the frac tion of warm 



molecular gas to the total amount of molecular gas. Using the CO survey of (jDame et al.ll200lh we integrate 
the amount of CO along the line of sight towards the SNR. We estimate the total column of molecular 
hydrogen using the canonical conversion N(H2)/Ico=2xl0 20 cm -2 (K km s -1 ) -1 . Though our apertures 
differ greatly, if the diffuse component is fairly uniformly distributed we estimate a warm gas fraction of 
0.005 and 0.15 towards SNRs G346.6-0.2 and G348.5-0.0 respectively. 



6. Conclusions 

The nature of six SNRs identified as having molecular colors in the GLIMPSE survey have been resolved 
spectrally by Spitzer IRS observations. Strong H2 S(0) to S(7) lines show clear evidence of shock interaction 
with dense gas. This doubles the number of known molecular remnants for which H 2 emission has been 
characterized. The excitation of H 2 is well fit by warm T~400 K and hot T~1500 K gas components. The 
OPR of the warm component deviates from equilibrium consistent with the passage of a slow shock through 
dense gas. 

Ionic fine structure lines are also present for all SNRs but have large variations in brightness between the 
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rcmnants. The presence of both strong molecular emission and ionic fine structure lines with large ionization 
potentials requires multiple shocks present within the observed aperture. In 3C 396 we see that the H2 
and ionic lines arc spatially separated, indicating changing physical conditions across the shock front which 
agrees with a multi-phase pre-shock medium composed of atomic, molecular and dense clump gas. 

Comparison of shock models shows that two C-shocks arc required to explain the observed H2 emission. 
Fitted shock parameters show one shock is of low velocity (~10 km s _1 ) into dense gas (^10 6 cm -3 ) while 
another shock of higher velocity (^50 km s _1 ) is propagating into lower density gas (~10 4 cm -3 ). Ionic lines 
require a fast J-shock into diffuse gas typically with a density ^10 3 cm -3 . We find physical conditions for 
these fast ionic shocks imply pressures similar to those found in the molecular C- shocks into dense clumps. 
For G349.7+0.2 we find extremely high ionic line pressures which are several times higher than for the H2 
emitting gas, possibly indicative of powerful ejecta shocks. 

Long slit spectroscopy of these SNRs isolates diffuse emission from that of the remnant. This is par- 
ticularly important for the H2 S(0) and S(l) lines as well as diffuse interstellar cooling lines such as [Si II] 
34.8/im. We find the diffuse H2 gas is relatively cool with a temperature of 100-300 K and a warm gas 
fraction of 0.5-15 % . 



We thank Achim Tappe for updating atomic data and fruitcful discussion on atomic Fe lines. This work 
is based on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion 
Laboratory under NASA contract. J.W.H. is grateful for the support of the IPAC staff and the Spitzer 
Visiting Graduate Student Program. 
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Table 1. Properties of Supernova Remnants 





Kes 69 


3C 396 


Kes 17 


G346.6-0.2 


G348. 5-0.0 


G349.7+0.2 


RA(J2000) 


18:33:01.89 


19:03:56.21 


13:05:32.75 


17:13:43.82 


17:15:04.90 


17:21:24.90 


DEC(J2000) 


-10:13:44.2 


+05:25:49.7 


-62:40:06.0 


-40:18:06.0 


-38:33:41.0 


-37:29:11.6 


diameter (') 


20 


8x6 


8 


8 


10 


2.5 


distance (kpc) 


5.2 


>7 


>9.7 


11 


13.7 


22 


physical size (pc) 


30 


>16xl2 


>23 


26 


39 


16 


velocity (km s _1 ) 


+80, +69 


? 


7 


-76 


22 


+ 16 


N H (10 22 cm" 2 ) 


2.4 


4.7 


3.6 


2.68 


2.8 


5.8 


r(9.7) 


0.95 


1.8 


1.4 


1.05 


1.1 


2.3 


A v (mag) 


13 


25 


19.3 


14.3 


15 


31 


References' 1 


1.2 


3,4 


5,11 


6,7,11 


5,10,11 


8,9 



"Re ference s (llWilsonl i 19721) ( 2~lYusef-Zadeh et alj d2003|), OlHarrus fc Slanj il999|), (4lQlbert et al.l 
j2003f), (51 ICaswell et al.l dl975f),(6lDubner et al I ||1993|) . ^Koraleskv et al.l J1998T) . (8~ lsiane et al.l toO^T 
(9) Lazcndic et al ■I i2005l) , (lOl lRevnoso fc Mangumf d2000ri (11) this work. 



Table 2. Observed Brightness of Emission Lines After Local Background Subtraction 











Intensity (ergs 


cm 2 s 1 sr 1 ) 








Transition 


A(^m) 


Kes 69 


3C 396 


Kes 17 


G346.6-0.2 


G348. 5-0.0 


G349.7+0.2 


H 2 S(0) 


28.219 


4.29(0. 28)e-6 


6.85(2.90)e-6 


4.54(0.28)e-6 


1.34(0.25)e-6 


3.26(0.27 


E-6 


1.54(0.15)e-5 


H 2 S(l) 


17.035 


4.14(0. 03)e-5 


3.15(0.41)e-5 


1.00(0.01)e-4 


4.26(0. 05)e-5 


1.27(0.04 


E-5 


8.18(0.10)e-5 


H 2 S(2) 


12.279 


4.51(0. 02)e-5 


4.02(0.41)e-5 


1.58(0.30)e-4 


4.60(0.07)e-5 


2.69(0.07 


E-5 


1.40(0.02)e-4 


H 2 S(3) 


9.665 


8.44(0.24)e-5 


2.56(0.19)e-5 


2.20(0.24)e-4 


1.32(0.02)e-4 


2.89(0.11 


E-5 


6.89(0.17)e-5 


H 2 S(4) 


8.025 


1.24(0. 16)e-4 


1.14(0.12)e-4 


4.56(0. 55)e-4 


1.20(0.05)e-4 


8.76(0.32 


E-4 


2.32(0.05)e-4 


H 2 S(5) 


6.910 


2.03(0. 10)e-4 


3.29(0.21)e-4 


1.28(0.09)e-3 


3.64(0.19)e-4 


2.99(0.04 


E-4 


4. 61(0. 07) e-4 


H 2 S(6) 


6.109 


1.06(0.05)e-4 


6.15(0.78)e-5 


3.80(0.32)e-4 


6.45(0.56)e-5 


8.87(0.41 


E-5 


l 


H 2 S(7) 


5.511 


3.15(0.09)e-4 


1.60(0. 14)e-4 


7.56(0.07)e-3 


2.22(0.05)e-4 


2.47(0.05 


E-4 


3.20(0.10)e-4 


[Fc II] 


5.35 


1.27(0.01)e-4 


2.82(0. 11)e-4 


1.55(0.07)e-4 




8.55(0.07 


E-4 


1.49(0.02)e-3 


[Ar II] 


6.98 










2.05(0.05 


E-4 


1.58(0.01)e-3 


[Ar III] 


9.00 














2.35(0.41)e-5 


[Ne II] 


12.8 


1.50(0.01)e-4 


1.10(0.08)e-4 


1.18(0.03)e-4 


3.20(0.08)e-5 


3.50(0.01 


E-4 


2.92(0.02)e-3 


[Ne III] 


15.5 


3.72(0.24)e-5 


1.51(0.08)e-4 


4.43(0.02)e-5 


2.47(0.28)e-6 


7.42(0.03 


E-5 


1.11(0.01)e-3 


[Fe II] 


17.9 


8.37(0.27)e-6 


3.55(0.04)e-5 


1.17(0.06)e-5 




5.88(0.05 


E-5 


3.77(0.02)e-4 


[S III] 


18.7 


4.19(0.26)e-6 


3.65(0.03)e-5 


8.83(0.34)e-6 


1.62(0.26)e-6 


1.59(0.03 


E-5 


2.43(0.03)e-4 


[Fe II] 


24.5 


3.44(0.23)e-6 


1.31(0.38)e-5 


2.76(0. 13)e-6 




1.62(0.03 


E-5 


1.90(0.03)e-4 


[S I] 


25.2 


5.64(0.29)e-6 




7.94(0. 25)e-6 


1.95(0.30)e-6 


1.77(0.02 


E-5 


1.17(0.02)e-4 


[Fe II] 2 


25.9 


3.84(0.02)e-5 


9.94(0.01)e-5 


4.54(0. 01)e-5 


7.00(0. 19)e-6 


1.80(0.01 


E-4 


9.82(0.05)e-4 


[S III] 


33.5 


1.52(0.04)e-5 


5.58(0.09)e-5 


1.94(0.03)e-5 


3.45(0.30)e-6 


3.35(0.03 


E-5 


4.02(0.03)e-4 


[Si II] 


34.8 


2.53(0.01)e-4 


2.12(0.05)e-4 


2.30(0.01)e-4 


6.93(0.05)e-5 


7.14(0.02 


E-4 


2.12(0.02)e-3 



1 Blended with strong PAH feature. 2 Blended with [O IV] 25.9/jm line. 



Table 3. De-reddened brightness of emission lines assuming an extinction curve of Ry = 3.1. 



Intensity (ergs cm 2 s 1 sr 1 ) 



Transition 


A(/im) 


Kes 69 


3C 396 


Kes 17 


G346.6-0.2 


G348.5-0.0 


G349.7+0.2 


H 2 S(0) 


28.219 


4.58(1. 51)e-6 


9.69(4. 10)e-6 


1.33(0.22)e-5 


1.49(0.30)e-6 


4.01(0.27 


E-6 


2.11(0.19)e-5 


H 2 S(l) 


17.035 


5.50(0.06)e-5 


5.70(0.75)e-5 


1.98(0.13)e-4 


6.01(0. 16)e-5 


1.81(0.44 


E-5 


1.88(0.39)e-4 


H 2 S(2) 


12.279 


5. 75(0.07) e-5 


8.03(0.81)e-5 


2.68(0.17)e-4 


6.83(0.11)e-5 


4.07(0.11 


E-5 


3.92(0.66)e-4 


H 2 S(3) 


9.665 


1.45(0.08)e-4 


1.47(0. 11)e-4 


9.09(0.50)e-4 


3.44(0.38)e-4 


8.18(0.30 


E-5 


5.40(0. 17)e-4 


H 2 S(4) 


8.025 


1.56(0. 14)e-4 


2.15(0.23)e-4 


7.59(0. 55)e-4 


1.75(0.07)e-4 


1.28(0.04 


E-4 


3.53(0.68)e-4 


H 2 S(5) 


6.910 


3.49(0.28)e-4 


4.61(0.30)e-4 


1.71(0.09)e-3 


3.58(0.86)e-4 


3.66(0.05 


E-4 


5.19(1.06)e-4 


H 2 S(6) 


6.109 


1.13(0.12)e-4 


9.36(1. 19)e-5 


5.10(0.32)e-4 


8.36(0.66)e-5 


1.14(0.05 


E-4 


l 


H 2 S(7) 


5.511 


3.44(0.21)e-4 


2.48(0.21)e-4 


1.00(0.07)e-3 


2.78(0.09)e-4 


3.21(0.06 


E-4 


5.34(0.62)e-4 


[Fc II] 


5.35 


1.61(0. 10)e-4 


4.47(0. 07)e-4 


2.22(0.10)e-4 




1.13(0.02 


E-3 


2.50(0.02)e-3 


[Ar II] 


6.98 










2.51(0.05 


E-4 


2.30(0.01)e-3 


[Ar III] 


9.00 














1.19(0.20)e-4 


[Ne II] 


12.8 


2.01(0.09)e-4 


1.96(0.15)e-4 


1.83(0.04)e-4 


4.44(0.09)e-5 


4.93(0.01 


E-4 


5.51(0.03)e-3 


[Ne III] 


15.5 


4.74(0. 04)e-5 


2.43(0. 13)e-4 


6.37(0.03)e-5 


3.24(0.36)e-6 


9.84(0.04 


E-5 


1.88(0.01)e-3 


[Fc II] 


17.9 


1.16(0.04)e-5 


6.71(0.06)e-5 


1.90(0.10)e-5 




8.59(0.08 


E-5 


7.62(0.02)e-4 


[S III] 


18.7 


5.79(0.36)e-6 


6.88(0.05)e-5 


1.44(0.06)e-5 


2.33(0.38)e-6 


2.32(0.04 


E-5 


4.90(0.03)e-4 


[Fe II] 


24.5 


4.23(0.29)e-6 


1.97(0.38)e-5 


3.76(0.17)e-6 




2.06(0.03 


E-5 


2.97(0.03)e-4 


[S I] 


25.2 


6.88(0.36)e-6 




1.07(0.03)e-5 


2.44(0.32)e-6 


2.23(0.03 


E-5 


1.80(0.02)e-4 


[Fe II] 2 


25.9 


4.66(0.03)e-5 


1.45(0.02)e-4 


6.06(0.02)e-5 


8.68(0.21)e-6 


2.25(0.02 


E-4 


1.49(0.01)e-3 


[S 111] 


33.5 


1.76(0.05)e-5 


7.42(0. 19)e-5 


2.41(0.03)e-5 


4.06(0.30)e-6 


3.97(0.04 


E-5 


5.50(0.04)e-4 


[Si II] 


34.8 


2.90(0.04)e-4 


2.76(0.07)e-4 


2.82(0.01)e-4 


8.06(0.05)e-5 


8.36(0.03 


E-4 


2.84(0.04)e-3 



1 Blended with strong PAH feature. 2 Blended with [O IV] 25.9/j.m line. 



Table 4. Comparison of IR Line Luminosities 



Luminosity (Lq) 


Kes 69 


3C 396 


Kes 17 


G346.6-0.2 


G348. 5-0.0 


G349.7+0.2 


H 2 S(0)-S(7) 


1.4 


2.5 


21 


7.2 


6.8 


68 


H 2 S(0)-S(2) 


0.5 


0.9 


6.3 


2.2 


1.5 


37 


H 2 S(5)-S(7) 


0.6 


0.9 


8.2 


2.3 


3.0 


13 


Ar 










1.4 


36 


Fe 


0.7 


3.8 


3.2 


0.4 


26 


448 


Ne 


0.6 


2.0 


2.3 


0.4 


9.5 


358 


S 


0.3 


2.4 


1.5 


0.3 


4.9 


201 


Si 


4.4 


0.5 


15 


5.4 


85 


704 



Table 5. Fitted Excitation Parameters to Observed Molecular Hydrogen Lines 





Kes 69 


3C 396 


Kes 17 


G346.6-0.2 


G348. 5-0.0 


G349.7+0.2 


N(H 2 ) roarm 


1.5e20 


3.0e20 


3.8e20 


4.7e19 


1.2e20 


2.8e20 


T wa rm (K) 


303 


265 


322 


528 


261 


467 


OPR„ arm 


1.3 


0.65 


1.4 


3.0 


0.4 


0.99 


N(H 2 ) hot 


5.0E18 


6.9E18 


2.6e19 


8.3e18 


2.4e18 


5.2e18 


T hot (K) 


1465 


1270 


1152 


1241 


1773 


1647 


OPRkot 


3.0 


3.0 


2.5 


3.0 


3.0 


3.0 
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Tablc 6. Combination of Two C-Shock Model Fits 





Kes 69 


3C 396 


Kcs 17 


G346.6-0.2 


G348. 5-0.0 


G349.7+0.2 


C-shock 1 














n (cm -3 ) 


10 5 


10" 


10" 


10" 


10' 


10" 


v s (km s -1 ) 


10 


10 


10 


15 


10 


10 


*s 


0.043(4) 


0.089(2) 


0.28(3) 


0.055(2) 


0.095(9) 


0.61(3) 


OPR. 


1 


2 


1 


3 


1 


1 


P s (dyn cm -2 ) 


2xl0~ 7 


2xl0 -6 


2xl0~ 6 


4xl0~" 


2xl0~ 7 


2xl0~ 6 


N(H 2 ) (cm" 2 ) 


4.3e20 


5.4e19 


2.0e20 


2.1e19 


1.0e20 


4.1E20 


C-Shock 2 














n (cm -3 ) 


10 4 


10' 


10 1 


10 5 


10 1 


10 s 


v s (km s~ 1 ) 


40 


30 


40 


50 


40 


40 


*s 


0.029(1) 


0.0091(1) 


0.094(1) 


0.004(1) 


0.024(1) 


0.0071(7) 


P s (dyn cm~ 2 ) 


3xl0~ 7 


1.5xl0~ 6 


3xl0~ 7 


4xl0~" 


3xl0~ 7 


3xl0~" 


N(H 2 ) (cm" 2 ) 


1.5e20 


2.6e18 


5.0e20 


3.5e24 


1.3e20 


7.1e19 



Note. — Errors to fitted shock models are found to be approximately 0.5 for gas 
density logio^, 5 km s — 1 for shock velocity v.g , and less than 10% for ^>.g. 



Table 7. Combination of C- and J-Shock Model Fits 





Kes 69 


3C 396 


Kcs 17 


G346.6-0.2 


G348. 5-0.0 


G349.7+0.2 


slow C-shock 














n (cm -3 ) 


10" 


10" 


10" 


10" 


10" 


10" 


v s ( km s~ 1 ) 


10 


15 


10 


15 


15 


10 


*s 


0.132 


0.046 


0.13 


0.053 


0.064 


0.54 


OPR 


2 


2 


2 


3 


1 


1 


P s (dyn cm -2 ) 


2xl0~" 


4xl0 -6 


2xl0~ 6 


4xl0 -6 


4xl0- 6 


2x10-" 


fast J-shock 
n (cm -3 ) 
v s ( km s - 1 ) 

*s 

P s (dyn era -2 ) 


10 3 
90 

0.104 
1.4xl0~ 7 


10 1 
130 
0.005 

3xl0~ 6 


10 3 

50 

0.19 

4xl0~ 8 


10 1 
150 
0.0053 
3xl0~" 


10 3 
60 

0.057 
6xl0~ 8 


10 4 
110 
0.10 

2xl0~ 7 



Table 8. Comparison of model fitting for single and multiple shocks 



Model 


One C-shock 


Two C-shocks 


Two n = 10 4 C-shocks 


C & J-shock 


SNR 


(no, Vs) [xl] 


(n , v s ),(n , Vs) [ Xr ] 


(no, v a ),(n , v s ) [ Xr ] 


(no, v 3 ),(n , Vs) [ Xr ] 


Kcs 69 


(10 5 ,20)[60] 


(10 5 ,10),(10 4 ,40)[12] 


(10 4 ,20),(10 4 ,50)[24] 


(10 6 ,10),(10 3 ,90)[20] 


3C 396 


(10 6 ,15)[56] 


(10 6 ,10),(10 5 ,30)[8.0] 


(10 4 ,20),(10 4 ,40)[23] 


(10 6 ,15),(10 4 ,130)[65] 


Kes 17 


(10 5 ,30)[23] 


(10 6 ,10),(10 4 ,40)[4.1] 


(10 4 ,20),(10 4 ,40)[5.4] 


(10 6 ,10),(10 3 ,50)[48] 


G346.6-0.2 


(10 6 ,15)[14] 


(10 6 ,15),(10 5 ,50)[7.0] 


(10 4 ,30),(10 4 ,50)[4.9] 


(10 6 ,15),(10 4 ,150)[0.9] 


G348. 5-0.0 


(10 6 ,20)[15] 


(10 5 ,10),(10 4 ,40)[2.4] 


(10 4 ,10),(10 4 ,40)[4.4] 


(10 6 ,10),(10 3 ,60)[12.4] 


G349.7+0.2 


(10 6 ,15)[16] 


(10 6 ,10),(10 4 ,50)[0.6] 


(10 4 ,20),(10 4 ,50)[6.3] 


(10 6 ,10),(10 4 ,110)[2.9] 



a The values in the brackets are reduced Xr- The best solution is in bold 
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Table 9. Physical properties derived from ionic line diagnostics 



(source) a 


[Fe II] ratios 


[Si II], [Fe II] 


[Nc II] lines 


[Nc III]/[Nc II] 


v s 


T-> b 

Pressure 


SNR name 


n„ (cm" 3 ) T (K) 


n (cm -3 ) 


v s (km s~ 1 ) 


v s (km s -1 ) 




(dyne cm -2 ) 


Kes 69 


multi-T or n e c 


500-1000 


150 


150 


150 


2.8xl0~ 7 


3C 396 


270 23000 


500-1000 


100-150 


350 


120,350 


1.8xl0~ 7 , 1.5xl0~ 6 


Kes 17 


100 7000 


500-1000 


150 


150 or 200 


150 


2.8xl0~ 7 


G346.6-0.2 




100-200 


70-100 


80 


80 


l.lxlO -8 


G348. 5-0.0 


multi-T or n e c 


500-1000 


160 


160 


160 


3.2xl0~ 7 


G349.7+0.2 


700 7000 


10 4 


250 


80 or 250 


250 


l.OxlO -5 



a Denotes which ionic lines were used to derive the density, temperature or shock velocity. b The pressure is estimated based on n D 
and V s derived from Si, Fe and Nc lines and from [Ne III]/[Nc II] ratio. c possible multi-tcmpcaturc and multi-density solution (see 
Fig. 8); the ratio of [Fe II] 17.9/5 fim yields a low density (< 400 cm -3 ) and low (<2000 K) temperature. 



Table 10. Fitted Physical Conditions of Diffuse H2 





3C 396 


G346.6-0.2 


G348. 5-0.0 


G349.7+0.2 


N ff2 (cm 2 ) 


2.1c20 


5.3el9 


4.3e21 


6.4e21 


T (K) 


155 


361 


118 


154 


OPR. 


= 3.0 


= 3.0 


= 3.0 


= 3.0 




1.2c20 


3.6cl9 


2.40e21 


3.6c21 


T (K) 


177 


509 


130 


176 


OPR. 


= 2.0 


= 2.0 


= 2.0 


= 2.0 




3.3cl9 




3.7c21 


1.5c21 


T (K) 


292 




222 


233 


OPR. 


0.65 




0.4 


0.99 



- 19 - 



a) Kes 69 





E 
4 



b) 3C 396 




PAH PAH 




c) Kes 17 



PAH PAH 



d) G346.6-0.2 



JJ 




PAH PAH 



Wavelength( J um) 



Fig. 1. — Observed IRS spectra after local background subtraction. Detected lines and PAH features are 
labeled. From top to bottom: (a) Kes 69, (b) 3C 396, (c) Kes 17, (d) G346.6-0.2, (e) G348. 5-0.0 and (f) 
G349.7+0.2. 
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H 2 S(7) S(6) S(5) 



Slow shock - 
Fast shock -: 




5.3|Jm X(nm) 8.7|jm 



Fig. 2. — Left: Spectral images for the short-low second order of IRS for 3C 396 between 5.1 and 8.7/im. 
Right: MIPS 24/xm image of SNR 3C 396 with the IRS SL2 slit overlaid. Arrow is drawn to indicate the 
direction along the slit. 




Fig. 3. — Left: Spectral image for the long-low first order of IRS for 3C 396 between 19.5 and 42.4/im. Note 
that the spectrum is severely degraded at wavelengths greater than 35/im toward the bottom of the spectral 
image. Right: MIPS 24/im image of SNR 3C 396 with the IRS LL1 slit overlaid. Arrow is drawn to indicate 
the direction along the slit. 
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G349.7 outside 



[Ne III] -> 

H 2 S(1) -> 

[Fe II] -> 
[S III] -> 




Fig. 4. — Left: Spectral image of the long-low second order of IRS for G349.7+0.2 between 14 and 21/iin. 
Diffuse PAH emission can be seen around the H 2 S(l) line. Right: MIPS 24^m image of SNR G349.7+0.2 
with the IRS LL2 slit overlaid. 



2000 



S(7) 



3 
2 

j=4 



/ 

/ 

S(6) / 
/ 

S(5) / 1-0 S(1) 
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S(4) 



/ 



S(3) / 
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v=0 v=1 
Vibrationol State 



Fig. 5. — Rovibrational energy level diagram for molecular hydrogen with rotational J states labeled. The 
S(0) to S(7) transitions available to Spitzer low-resolution IRS observations are indicated with solid lines. 
All 'S' transitions are from the upper to lower J state with AJ=2. The v=l-0 S(l) transition commonly 
observed at 2.12^m is plotted as a dashed line. Note that the S(7) transition traces a higher upper excitation 
state than the 1-0 S(0) and S(l) lines commonly observed in the NIR. 
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Fig. 6. — Boltzmann diagram for the excitation of H2. From top to bottom: (a) Kes 69, (b) 3C 396, (c) 
Kes 17, (d) G346.6-0.2, (e) G348. 5-0.0 and (f) G349.7+0.2. Triangles mark the column density of the upper 
J state obtained from the de-reddened brightness with error bars. The best-fit two-component model is 
ovcrplotted as a solid dark line. A dashed line shows the contribution of the warm component; a dotted line 
shows the contribution of the hotter component. Parameters for the fits are given in Table [5] 
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01 2345678 
Energy of upper level (1000 K) 



Fig. 7. — Two C-shocks arc fitted to the observed excitation of H2. From top to bottom: (a) Kes 69, 
(b) 3C 396, (c) Kes 17, (d) G346.6-0.2, (e) G348. 5-0.0 and (f) G349.7+0.2. The total of two fitted shock 
components is overplotted as a solid line A dashed/dotted line shows the contribution of the slower/faster 
shock. Parameters for the fits are given in Table [6l 
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2 4 6 G 02468 02468 

Energy of upper level (1000 K) Energy of upper level (1000 K) Energy of upper level (1000 K) 



Fig. 8. — A comparison of fitted shock models to excitation of H2 by a combination of shocks for Kes 69 
[left], G346.6-0.2 [center] and G349.7+0.2 [right]. From top to bottom: a single C-shock, two C-shocks with 
different densities and velocities, two C-shocks with the same density of 10 4 cm' 3 , a combination of one 
C-shock and one J-shock. Table [8] has the fitted values and reduced chi-squared values. In cases of multiple 
shocks the slower shock is plotted with a dashed line and the faster shock with a dotted line. The total 
contribution of the two shock components is plotted as a solid line. 
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Fig. 9. — Contour plots of of [Fc II] diagnostic line ratios using 5.3, 17.9 and 26^m lines. The observed ratios 
of G349.7+0.2 (square), 3C 396 (triangle) and Kes 17 are marked. The two ratios for Kes 69 (two thick 
dotted lines) and G348. 5-0.0 (two thick solid lines) do not converge to the same temperature and density 
ranges which may be due to uncertainties of collisional strengths at low density and temperature or presence 
of multi-temperature components. 
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Vs (km s ) 

Fig. 10. — Predictions from the models of Hartigan et al. (1987) of the [Ne III] 15.5/[Ne II] 12.8yum ratio as 
a function of shock velocity for densities of 100 and 1000 cm~ 3 . The observed values for each remnant are 
marked as dashed lines. 
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Fig. 11. — Boltzmann diagram for the excitation of diffuse H2 observed toward SNRs in our sample. The 
SNRs are color-coded with blue for G346.6-0.2, red for 3C 396, purple for G349. 7+0.2, and yellow for 
G348. 5-0.0. 



